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ABSTRACT 

The first spherical accretion model was developed 55 years ago, but the the- 
ory is yet far from being complete. The real accretion flow was found to be 
time-dependent and turbulent. This paper presents the minimal MHD spherical 
accretion model that separately deals with turbulence. Treatment of turbulence is 
based on simulations of several regimes of collisional MHD. The effects of freezing- 
in amplification, dissipation, dynamo action, isotropization, and constant mag- 
netic helicity are self-consistently included. The assumptions of equipartition 
and magnetic field isotropy are released. Correct dynamics of magnetized flow 
is calculated. Diffusion, convection, and radiation are not accounted for. Two 
different types of Radiatively Inefficient accretion flows are found: a transonic 
non- rotating flow (I), a flow with effective transport of angular momentum out- 
ward (II). Non-rotating flow has an accretion rate several times smaller than 
Bondi rate, because turbulence inhibits accretion. Flow with angular momen- 
tum transport has accretion rate about 10-100 times smaller than Bondi rate. 
The effects of highly helical turbulence, states of outer magnetization, and dif- 
ferent equations of state are discussed. The flows were found to be convectively 
stable on average, despite gas entropy increases inward. The proposed model has 
a small number of free parameters and the following attractive property. Inner 
density in the non-rotating magnetized flow was found to be several times lower 
than density in a non-magnetized accretion. Still several times lower density is 
required to explain the observed low IR luminosity and low Faraday rotation 
measure of accretion onto Sgr A*. 

Subject headings: accretion, accretion disks — MHD — turbulence - Galaxy: 
center 
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INTRODUCTION 



Dynamics of magnetized accretion flows is a major topic of astrophysical research. The 
problem can be solved with two different approaches: numerical and analytical. Each of 
them has specific difficulties, so these methods can be applied together for a better result. 

Realistic num erical simulatio ns require a lot of computational time to model even 
the isotropic case (ILazarianl 120061 ) . Convergence of properties of the isotropic turbulence 
i s reached only when computational doma in has more than f024 cells in each dimension 
( jLadeinde fc Gaitondell2004l ; iBiskampl 120031 ). Non- isotropic simulations of this size were not 
performed. It is also very difficult to model the system with large range of scales. The 
system then possesses va stly different tim escales. Existing simulations of accretion flows 



object (IHawley &: Balbus 



2002 



are eit her axisymmetric ( McKinnevi 120061) or consi der a rather small domain close to the 



I gumenshchevi 120061 ) . In addition, simulations should be run 



for sufficiently long time or several runs should be made to obtain average quantities, e.g. 
accretion rate, power of emitted radiation. 

Analytical models do not suffer from a need to average, if they are based on averaged 
quantities. However, to build a reasonable model is itself difficult. No unified method exists 
to combine insights in physics and mathematics into a perfect analytical model. That is why 
the zoo of approximations of astrophysical flows is so huge. 



In particular. 



( jNaravan &: Yi 



many 



metric treatme nt (Bondi 1952 ; Meszarosj 



analytical treatments were devised for accretion: spherically sym - 



9751 : ICoker fc Meliall2000l : iBeskin fc Karpovll2005h , 
standard disk (IShakura &: Sunyaevi 1 19731 ). Advection-Dominated Accre tion Flow (ADAF) 



19951 ) with its variation H ot Luminous Accretion Flow (jYuanl 120011 ). Adia- 



batic Inflow- Outflow Solutions (ADIOS) (IBlandford fc Begelmanlll9991). Convection Domi- 

nated Accretion Fl ow fCDAF) (jNaravan. Igumenshchev. fc Abramowiczll2000l : lOuataert fc Gruzinov 
2000), Jet-ADAF (lYuan. Markoff. k Falckel l2002h . They are aimed to describe essentially 
the same process: axisymmetric plasma inflow onto a compact source. Some models include 
the effects the others miss. Energy transport in CDAF, outflows in ADIOS are the examples. 
Some effects are not treated properly in any approximation. 

Magnetic field is a main source of uncertainty and mistakes in theory of accretion flows. 
Two assumptions are usually posed to incorporate i t into the model. Firs tly, magnetic field 
is considered to be isotropic flCoker fc Melial 12000; iNaravan &: Yil 119951) . Then magnetic 
pressure and magnetic energy density may be put JNaravan fc Yilll995h into the dynamical 
equations. Secondly, the ratio of magnetic field energy density to gas thermal energy density 
is set to constant. This is called thermal equipartition assumption. These two ideas are at 
least unproven or may even not work. Magnetic field is predominantly radial in spherical 
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( IHawley fc Balbus 



i nflow (|Shvartsmanll97ll) because of freezing- in condition and predominantly toroidal in disk 



200 



1) because of magnetorotational instability. 



In a good model direction and strength of the magnetic field should be determined self- 
consistently. Non-isotropy of magnetic field re quires special dynam ics. Dynamical equations 
were partially derived more than 20 years ago (Sc harlemann 1983), bu t did not receive much 
attention or were even considered erroneous (iBeskin fc Karpovil2005l ). 



Such a model may offer a natural explanation of certain accretion patterns. Accretion 
onto Sgr A* gives an excellent opportunity for testing. Ou r Galaxy is proven to host a 
Supermassive Black Hole (SMBH) named Sgr A* in its center (IGhez et al.ll2003l ; IShenll2006l ). 
This black hole accretes m atter and emits radiation with characteristic low-luminosity spec- 
trum ( Narayan et al.l[l998l ) . This spe ctrum was satisfactory explained w ith the combination 
of two models: jet or non-thermal ([Yuan. Quataert. fc Narayan! 120031 ) radio-emission and 
X-Rays with IR radiation coming from conventional ADAF flow. However, the large number 
of free parameters allows one to fit any spectrum well. Model with no free parameters left 
is an ultimate goal of the ongoing study. 

Partial progress in building a self-consistent accretion model is made in this paper, 
which is organized as follows. Averaged spherical MHD model with turbulence is devised in 
Section [2j Approximate model employs the characteristic length scale about the size of the 
region of interest. Coefficients are taken from several hydrodynamic and MHD simulations. 
External sources sustain turbulence at large radii, whereas turbulence is self-sustained in the 
converging flow at small radii. Necessary boundary conditions are discussed in Section [3] for 
general flow and for Sgr A*. Results in Section H] are followed by the discussion of the model 
in Section [51 Observational implications in Section [6] are supplemented with prospects for 
future work and Conclusion in Section [7J Paper has several appendices. 



2. SPHERICAL MODEL 



I base all calculations on Magneto Hydrodynamic system of equations (ILandau. Lifshitz fc Pitaevskii 



19841 ) . The viscous terms are retained where they do not vanish in the limit of vanishing vis- 
cosity. The quantities in the following equations are fully dependent on time and coordinates. 
General mass flux equation reads 

^ + V( P V) = 0, (1) 

where V is fluid velocity. Force balance is described by Navier-Stokes equation 
<9V , Vp [B x [V x Bll 

_ + ( VV)V^-V^- ' I " + ,AV, (2) 
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where g is gravitational potential, v is kinematic vis cosity. The last term is r esponsible for 
finite energy dissipation through Kolmogorov cascade (ILandau &: Lifshitzlll9871 ). Momentum 
equation is a combination of equations ([TJ) and (J2J) 



<9 



<9t (9x fe 
Energy equation 

d ( oV 2 B 2 \ / /V 2 \ 1 



+ i/(AV), 



(3) 



at 



8tt 



47T 



VISCOUS 



J (4) 



includes information about the equation of state. Here e is gas internal energy density, 
w = e + J dp/ p is gas specific enthalpy. Viscous term is responsible for diffusion. Magnetic 
field evolution is described by induction equation 



<9B 

— = V x [V x B] + u M AB 



(5) 



with magnetic diffusivity vm- Magnetic field is solenoidal as well as incompressible random 
velocity field: 

VB = 0, Vu = 0. (6) 



2.1. Dynamics 

Spherical accretion is the simplest pattern of all symmetric setups. We need to solve 
the basic model first to move then to a more realistic pattern. Construction of the minimal 
maximally symmetric model is the subject of the following study. 

I employ the natural for the problem spherical coordinates (r, 6, <fi) and average over 
angular variables (0,<p). The results depend only on the radial variable r and not on time 
t in the assumption that angular averaging is the same as time averaging. I need now to 
determine the essential quantities and derive the closed system of equations on them. 



Essential quantities of a non-magnetized solution in iBondil ( 119521 ) are the inflow speed 
v (r), density p(r), and temperature T(r). Turbulent magnetized case requires several more. 
As I release the assumption of isotropy, there are two special directions: along the radial vec- 
tor e r and perpendicular to the radial vector. To describe realistic Magneto Hydrodynamic 
turbulence, I need at least 6 quantities: squares of radial and perpendicular magnetic fields 
Bl and B\ , squares of radial and perpendicular random fluid speeds u 2 and u 2 ± , characteris- 
tic length scale L, and dimensionless magnetic helicity £. The last quantity will be described 
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in detail in the corresponding subsection 12.41 For simplicity I consider random velocity to 
be isotropic and denote it as u(r). 

Total velocity of a fluid parcel 

V(r,e,</>,t) = v(r)er + u(r,e,<l>,t) (7) 

is a sum of averaged inflow speed v(r) and instantaneous random velocity u(r, 9, <p, t), where 
by definition angular average of turbulent velocity vanishes 

J u{r,9,(f),t)dn = 0. (8) 

General continuity equation (JTJ can be averaged with the aid of equations (jUJ) and (jSJ) to 

Anp(r)v(r)r 2 = M, (9) 

where M is the mass accretion rate. 

I derive the averaged force equation from general momentum equation d3J). Tensor 
pViVk averages out into the diagonal form pv 2 5 rr + pu 2 5ik/3- Because there are no sources 
of magnetic field (eq. [fi]) and sphe rical geometry is assumed, no regular magnetic field 



exists. Following IScharlemannl (119831 ). I add -B r VB/(47rp) to the radial magnetic force F r = 
[B x [V x B]] r /(47rp), average over the solid angle, and then set = B± and Be = B±. 
Cross-terms with (BqB t ), (B^B r ), and (B^Bg) cancel on average over the solid angle. Finally. 
I obtain 

P _ (r A B 2 Y r (r 2 Bl)' r 



, > ( 10 ) 
7rpr 4 Airpr z 



omit bulk viscosity term that 



for the magnetic force. I denote by ()', radial derivatives. 

results from uAV. Paczynski-Wiita gravitational potential fjPaczhynski fc Wiitalll980T ) 

(id 



2(r - r g ) 

is used to imitate the effects of General Relativity, where 

2GM 



c 2 



(12) 



is a Schwarzschild radius of an object with mass M. I take gas pressure to be that of an 
ideal gas p = pRT/p, where p is a mean molecular weight. Combining all the terms, I come 
to the averaged force equation 

, , V 2 R (pT)> (pu 2 )> (r 2 Bl)> (r*B?)' r 
2{r — r g ) z /i p dp 4:7rpr z 87rpr 4 
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Averaged energy advection equation can be derived directly from general energy equa- 
tion (j3J). Enthalpy term should include contribution from random fluid motions as well as 
from gas. Isotropic random motions of fluid exert isotropic pressure p ran d = pu 2 /3 and have 
the internal energy density £ ran d = u 2 /2. Total enthalpy w is 

, i?r(/ e q e (T) + M(r) + l) 5 2 

w = w gas + w rand , where w gas = and w rand = -u . (14) 

fi 6 

Fractions of electrons f e ~ 0.54 and ions f\ ~ 0.46 are calculated for a gas with twice 
solar abundance of el ements. Such high concentration of helium and metals was assumed by 



Baganoff et al.l (120031 ) for spectrum fitting of Sgr A*. Correspondent mean molecular weight 



is \x ~ 0.7g cm 3 . Integral heat capacity per par ticle a P {T) and aAT) are different for electrons 



and ions. Ions are n on-relativistic down to r g (INarayan &: Yilll995l ). Therefore Oj(T) = 3/2 



General expression (jChandrasekharl 119571 ) should be used for thermal relativistic electrons 
a e (T) = e- 1 (3fr 3 (e- 1 ) + ir 1 (e- 1 ))/(4 J RT 2 (e- 1 ) - l). Here 6 = kT/m e c 2 is dimensionless 
temperature, K X (Y) are modified Bessel functions of the second kind. Expression for non- 
relativistic enthalpy is 

5RT 5 2 

wnr = — h -u . (15) 

Zfi 6 

It is valid in the limit 9< 1. Time derivatives in energy equation (jl]) vanish under averaging. 
Equation takes the form Vq = 0, where q is the energy flux. Part of flux proportional to 
random velocity u averages out, because turbulence is incompressible and u is zero on average 
(eq. [8]). Applying continuity relation (jHD, I finally obtain 

VV r + WT^ v?+ w r + 7^{ — ) =0, (16) 



2(r - r g ) 2 ' ~ T ' 2n \ p 

where again By = B 2 = B\. I assumed the term J[Bx [u x B]]df2 to also be zero along 
with all viscous energy transfer terms. I limit this study to Advection Dominated flows by 
deliberately cutting off diffusion and convection (see Appendix [Cj . 

Subtracting force equation ( |T3l) from energy advection equation (flEI) I get the heat 
balance equation that reads in non-relativistic limit 



similar to entropy conservation in hydrodynamics. Work done by gas is represented by 
— p' r j pT. The first term has exactly the form of the second, if I make the substitution of the 
mean square particles velocity 



vl = (1* 



Work done by the magnetic field enters the expression as derivatives of p and r in the 
magnetic part. 
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2.2. Evolution of Turbulence 



Dynamics is the only part of ideal Bondi problem (lBondilll952l ). In reality, flow always 
has some small scale turbulence that exerts back-reaction on the mean flow. The magnitude 
of back-reaction terms should be determined from additional equations that describe the 
evolution of random magnetic field and fluid motions. Since no complete theory of turbulence 
exists, I make a lot of approximations. The model is adjusted to agree wit h the results o f 



several numerical simulations. I also apply analytical tests similar to that in lOgifvid (120031 ) 
to assure the model reproduces the basic properties of observed turbulence. 

I need non-ideal induction equation (j5J) and Navier-Stokes equation (j2J) to derive how 
turbulence evolves. My goal is to compound reasonable equations on average squares of 
radial magnetic field Bf, perpendicular magnetic field Bj_, isotropic velocity u 2 . I also need 
equations on characteristic length scale of turbulence L and dimensionless magnetic helicity 

£■ 

Radial part of induction equation <^ easily gives the equation on Bf, when the former 
is multiplied by 2B r and averaged over the solid angle: 

OB 

2B r-^r- = 2B r [Vx [ve r xB}} r + 2B r [V x [u x B]] r + 2u M B r (AB) n (19) 

where indices () r without primes denote the radial parts. The left-hand side vanishes as 
all time derivatives. The first term on the right-hand side represents the uniform increase 
of magnetic field due to flux freezing. I combine it with the continuity equation ([9]) to 
eliminate v derivatives. The second term is the dynamo action. It cannot be easily averaged. 
Characteristic turbulence length scale L may be used to approximate derivatives 

dBi Bi dui Ui . . . 

a — ~ -r^k and - — ~ — e k , e k - unit vector. (20) 
oxk L oxk L 

Then we arrive at dynamo action with characteristic timescale Td yn = CB U T~eda about eddy 
turn-over time r c dd — u/L. The averaged expression is quadratic in magnetic field. I take 
coefficient to be cb u \ at any Bf and cb U 2 at any BiB^ with i ^ k. The final form of the 
dynamo term reads 2B r [V x [u x B]] r = (csuiBf + CBuiB r {Bg + B^))u/ L, and characteristic 

B r = v^Bf and Bq = B(j> = B± = (21) 

should be taken. The last term on the right-hand side of equation (1191) represents magnetic 
field dissipation. Dissipation term %AB of inducti on equation (Q is macroscopic in tur- 
bulence even for vanishing magnetic diffusivity vm ( Biskamp 20031 ) . I approximate radial 



dissipation to have a timescale Td yS s = cbb^At about Alfven timescale ta t = VAr/L. The av- 
eraged expression is also quadratic in magnetic field. I take coefficient to be cbbi at any Bf 
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and c BB2 at any BiB k with % ^ k. Finally, v M B r (AB) r = va{c BB \B 2 + c BB2 (B e + B^)B r )/L. 
Collecting all the terms, I obtain 

v d(B 2 r 4 ) = -(c Bul B 2 + 2c Bu2 B r B ± )u + (c BB1 B 2 + 2c BB2 B r B ± )v Ar 
r 4 dr L 

for the radial magnetic field in the absence of external energy sources. 

Perpendicular part of induction equation for example 8 part, gives the equation 
on Bq when equation (jSJ) is multiplied by Bg and averaged over the solid angle. The flux 
freezing condition for perpendicular field is different from that for radial field: Bgvr = const 
represents perpendicular flux freezing. I repeat the calculations made for radial field B r to 
find dynamo and dissipation terms. Dynamo term takes form (cBuiBg+CB U 2Be(B ( j > +B r ))u/L. 
Dissipation term is v Ae(c BB iB 2 + c B B2(B<t> + B r )B g )/L with perpendicular Alfven timescale 
for dissipation. Here I take B\ = Bg = BgB^ = Finally, I obtain 

2 2 9 ( B l \ -((CBUI + CBu2)B]_ + C B u2B±B r )u + ((c BB1 + c B m)B\ + c BB 2B±B t )va± 

dr \p 2 r 2 J L 

(23) 

where continuity equation §§§ is used. Radial Va t and perpendicular va± Alfven speeds and 
random velocity u are 

V Ar = t= , v A ±_ = 7= - , u = Vv?. (24) 
V 47r P V 47r P 

Coefficients c Bu i, cb U 2, cbbx, cbb2 are yet to be determined. 

Evolution equation for squared random fluid velocity u 2 can be found from momentum 
equation ([3]), when it is multiplied by 2u and averaged over the solid angle. Potential energy 
and pressure terms average out and only three terms are left 

(, x V(pV)\ u[B x [V x Bll , . 

2u V V V + KH ' = 2^ [ - ^ + 2uz/Au. 25 

V P J 4vrp 

I apply the same averaging procedure as for magnetic field evolution equations (1221) and (I23p . 
The final result is 

2/3 9 ( \ - CuuVl3 ~ ( C uBlVA + (2CuBl + C uB2 )v 2 A± + 2c uB2 (v A V 'A±))u 

VP dr Vp 2 / 3 )~ L ' 1 ] 

with additional three coefficients c uu , c uB \ and c uB2 . Some of these and other c xx -\ike coeffi- 
cients can be taken from numerical simulations of isotropic turbulence, some of the m can be 
i nferre d from analytical tests. They may not simply be set to convenient values like lOgilvie 
f l2003h did. 
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2.3. Correspondence to Numerical Simulations 

Isotropic turbulence is studied quite thoroug hly in numerical simulations. Some results 



are reproduced by a number of researchers (see iBiskarnpl (120031 ) for the review). That is 
why we may believe in these results and base a model on them. Three simulations of 
different turbulence regimes can provide four conditions that let us uniquely determine four 
combinations of coefficients c xx . These regimes are decaying HD turbulence, decaying MHD 
turbulence, and dynamo growth of small seed magnetic field. I assume then that c xx are 
constants independent of regime and extend the derived model to any anisotropic case. 

Let me consider my model in isotropic incompressible case of box turbulence. In these 
settings B 2 = Bg = B\. Squared magnetic field B 2 equals B 2 = 3B 2 . Transition to the 
co-moving frame of averaged inflow in turbulence evolution equations (1221) . ff23l) . (126]) is done 
by stating d/dt = —vd/dr. Now I should write time derivatives instead of radius derivatives 
and set r = const, since matter is not moving anywhere from the box. I obtain equations of 
evolution of isotropic turbulent Alfven speed va and isotropic turbulent velocity u: 

/ 2 \' _ CuBv\ U ~ CuuU 3 / 2 y _ CBuV A u ~ C B Bv\ / ov n 

\ U It ~ T ' \ V A) t - T • W 



Here va = v B 2 /y/Airp and p = const. Coefficients with hats are 

, o - CBB1 + 2C 5£2 /r>Q\ 

CBu — CBul + *CBu2, c BB — ~ ^= , l^oj 

Cuu Cuul i CuB CuBl ~r C U B2 

in terms of previously defined c xx . 

I have a freedom to set L, because it enters the equations only in combinations c xx /L, 
but c xx are not yet determined. For simplicity of further derivation I take L(r) to be the 
effective size of energy containing eddies for isotropic incompressible turbulence: 

POO POO 

2 _ / L. |2 .7 „„ , 2 _ / I |2 



u = \uk\ dk and v A = \vAk\ dk. (29) 

J2tt/L J2tt/L 

Isotropic decay of hydrodynamic turbulence is the simplest simulation. The convenient 
constant of decay is Kolmogorov constant Chd- It is defined as 

C HD = E k e'h- 2 ^ with e = -±(j\ and E k = l -^f, (30) 
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where Ep, is energy spectrum, e is a dec ay rate. Kolmogo rov constant was found to be 



Chd ~ 1-65 in the large set of simulations (jSreenivasanl 1 1 9951 ) . I substitute this number into 



equation (13"01 and evaluate the first integral in equation (}2"9"]) to find 

47T 



(3^)^ 

for isotropic equations (12"7|) . 



1.14 (31) 



Isotropic decay of magneto hydrodynamic turbulence gives two conditions. MHD Kol- 
mogorov constant is defined similarly to HD case equation (130]) as 

C MHD = E k k?<*e-W with e = ~(^±^j and E k = ^ ( 32 ) 
MHD turbulenc e is more diffic ult to model numerically, but the value of Cmhd ~ 2.2 is 



rather rigorous (IBiskampI 120031 ) . In addition, kinetic energy was found to decay in exactly 
the same rate as magnetic energy. Evaluation of the sum of two integrals (129]) with definitions 
( B2J) and known Cmhd yields 

/ 2 \ 3/2 

cbb - c Bu = Cuu - c uB ~ 27T — w 1.05. (33) 

V^MifO/ 

Dynamo simulations explore the regime v\ <C u 2 . Exponential growth of small magnetic 
field corresponds to some value of coefficient cb u in equations (l27j) as 

5 2 oc exp f . (34) 

External driving is purely mechanical for v\ <C u 2 , so external source of magnetic field does 
not alter the picture of field amplification by dynamo. Characteristic length scale in dynamo 
simulations is usually the size of energy containing eddies L consis tent with definition (1291) . so 



renor malization of length scale is not required. Older simulations (IKida. Yanase. fc Mizushima 



1991 ) have found b = 0.39 that corresponds to &b u ~ 0.61. Later results (ISchekochihin et al. 



20041 ) indicate a bit higher value &b u ~ 0.7 that I will use for my model. Finally, 



cbu = 0.70, c B b = 1-75, (35) 
c uu = 1-14, c u b = 0.09. 

The values of four c xx (eq. [35]) are not enough to obtain all seven coefficients c xx in 
equations f[2"2"]) . fl2"3"]) . fl2"Bl) with definitions (T2"g]) . However, the application of common sense 



- 11 - 



analytical conditions to non-isotropic system of equations puts some additional constrains 
on c xx that allows me to complete the model with as little guessing as possible. 

Analytic tests are described in Appendix |A] This completes the derivation and verifi- 
cation of turbulence evolution equations (122]) . (1231 . (|26|) with coefficients 

c BBl = 3.03, cbb2 = 0.00, c Bul = 0.41, c Bu2 = 0.29, (36) 
c uu = 1.14, c uB1 = 0.09, c uB2 = 0.00 

that I obtain summarizing equations (B5jl. (j55jl . flATTj) . and flAT2|) . However, not all major 
effect have been included so far. 



2.4. Magnetic Helicity 



Certain correlation called "magnetic helicity" may strongly influence magnetic field 
dissipation. This quantity is defined as 



H 



(A B)dV, 



(37) 



where A is a vector potential with a defined gauge condition (jBiskampll2000i ) . Time derivative 
of magnetic helicity is very small compar ed to the time derivative of magnetic energy in high 
Reynolds number astrophysical plasma (jBiskampl 120031 ) : 



dH E 



M 



dE]\j H 



<C 1. 



(38) 



Constancy of magnetic helicity defines the rules of selective decay. Magnetic energy 
Em decays in free turbulence down to non-zero value, allowed by constant magnetic helicity 
H = const. The final force-free configuration has zero random k inetic energy Er and has 
aligned current density and magnetic field j ff B ( jBiskampl 120031 ) . 



However, the derived system of turbulence evolution equations ( 1A1I) and, therefore, 
equations (j22j) . (1231 . (|26l) cannot handle selective decay. Decay of magnetic energy must 
be modified in order to have the transition to zero dissipation rate at certain va t and vai_ 
as a function of magnetic helicity H. First, I should employ the proper magnetic helicity 
constancy. Then I should quantify the relation between critical VAr, va±, and H. 

Let me consider the region S that evolves together with the mean flow of fluid. This 
region has the constant angle boundaries 9 = const and <fi = const. Its radial elongation L r 
scales as inflow velocity: L r oc v . The region S contains constant mass m = const of matter, 
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because matter flux through its boundaries is zero by definition. If I neglect diffusion by 
random velocity, frozen magnetic field lines do not move through t he boundaries of the 



region. Because of this, magnetic helicity in S is constant H = const (IBiskampI 120031 ) . 
The simplest order of magnitude relation between magnetic energy Em and H is 

E M L H = H = const (39) 



in the region S, where Lh is magnetic helicity characteristic length scale (IBiskampI 120031 ) . 



As magnetic field decays in turbulence, Lh grows according to equation (1391) . 
I can parametrize Lh to be a fraction of L : 

L H = £L. (40) 

Volume of the region of interest S is 

V = 

P 



V=™ (41) 



with m = const. Total magnetic energy Em is 

V 



E M = —(B 2 r + 2Bi). (42) 



■57T 



I substitute relations (j4"0|) . fj4TT) . and f|42|) into equation (1391) and use the definitions (124|) of 
Alfven velocities to come to 

L(v\ r + 2v\ L )i = const. (43) 

Now I need to include £ into the turbulence evolution equations (122]) . (j2"3"j) . fl26|) so that 
they can handle selective decay. The natural limit of Lh growth is the characteristic size 
of energy containing eddies L. So regime ^ < 1 corresponds to non-helical turbulence and 
regime £ ~ 1 to turbulence, where magnetic helicity significantly inhibits dissipation. Regime 
£ ^> 1 does not occur. The basic way to modify the equations is to decrease by a smooth 
multiplier / (£) < 1 magnetic field decay rate. For qualitative agreement with experiment 



(IBiskampI 120031 ) I can employ 

/(£) = exp(-£), (44) 

what means that magnetic energy dissipation timescale becomes exp(£) times larger. Terms 
with both u and one of VAr and va±_ in magnetic field evolution equations (!22j) . (|23|) do 
not need to be modified, since random velocity energy decays to zero and these terms do 
not matter. However, I multiply the term with both random velocity and Alfven speed in 
turbulent velocity evolution equation (f2"oT) by exp(— £) to make random velocity u decay to 
zero. 
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2.5. System of Equations with Source Terms 

With only minor corrections, the final system of equations can be written down. In 
general, turbulence has external sources of energy that sustain finite magnetic and kinetic 
energies even in case of box turbulence. I can add source terms to incompressible system 
OAip and consequently to the system of compressible equations 



System AA1H with coefficients (I35I) and (1361) . modifier (jHJ), and source terms reads 
d{v\ r ) _ (0.70v Ar + 0.58(v A± -v Ar )v Ar )u- 3.03v Ar exp(-£) v z p 

~ir- 1 + Co T' (45a) 

d(v 2 A± ) _ (0.70v 2 A± + 0.29(v Ar -v A± )v A± )u-3.03v A± exp(-O v 3 p 

L + Cl L' (45b) 

d(u 2 ) _ 0.09(v 2 Ar + 2v 2 A± )uex V (-Q - l.Uu 3 v} 
dt L +C2 L' { ' 

where v p is the mean square particles speed (eq. [18]) and c , Ci, and c 2 are dimensionless 
coefficients. These coefficients determine the rates of external energy input into turbulent 
fields. 

I denote by a the ratio of total turbulent energy to thermal energy: 

Ek + Em , 3i?T v 2 p u 2 v\ 2 

a = — , so that a— = = — + — + v A± . (46) 

Unlike conventional plasma magnetization, magnetization a with definition ( 1461) includes the 
energy of random fluid motions. 

In the dynamic equilibrium of constant VAr, va±, u and known £ system (j4"5j) gives three 
algebraic equations for ratios VAr/v p , va±/v p , and u/v p as functions of Co, ci, and c-i- To 
estimate cq, c\, and c 2 I take stationary driven isotropic turbulence with kinetic energy Ex 
equal to magnetic energy E M . Isotropic turbulence of interest has v Ar = va± = u/y/3. Such 
turbulence occurs far from the central object, where outer magnetization is a constant a^. 
Solving system (Jl5|) I obtain using equation PBT) 

Co 



ci w 0.124a^ /2 , c 2 = 3c w 0.371a^ /2 (47) 



in case £ = 0. I apply these values even to turbulence with £ > 0. Total external energy 
input Q + into Ek and Em is 



Q + « 0.742a^. (48) 

This energy adds up to thermal gas energy after being processed through turbulence. How- 
ever, I do not adjust my dynamical equations f|T3|) and ffTBT) for Q + . I self-consistently omit 
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external heating and radiative or diffusive cooling. This omission is physically justified suf- 
ficiently far from the central object, where cooling Q_ balances external heating Q + . It is 
also justified in the inner region, where both Q + and Q_ are negligible compared to the 
internal driving and energy advection. Internal driving represents bui ld-up of self-sustaine d 



turbulence in a converging flow due to conservation of magnetic flux (ICoker fc Melial 120001 ) . 



Only the size L of energy containing eddies should be specified to complete the derivation 
of closed system of equations. In the case when energy input Q + does not matter, the problem 
has only one relevant scale that is the size of the system r. Therefore, I can set L to be the 
fraction of radius 

L = 7r (49) 

with the proportionality constant 7 about unity. However, energy input from external sources 
Q+ is relatively large far from the central source. This causes medium with constant Q + , 
constant v p , and constant to have constant size of largest eddies 

L = Loo = const (50) 

because of equation fj4*8|) . This equality holds for radii larger than some tq « Loo/j. I 
introduce a function with a smooth transition from relation f|4"9]) for r r to relation (1501) 
for r ^> r : 

L(r)=L OQ (l-exp(-^-)) (51) 



This completes derivation and verification of 8 equations (jUJ), (ITB1) . (fTB"j) . f[2"2l . (1251) . f[2"B"j) . 
( 143|) . (l5TT) with coefficients (|35|) . (!36|) . and (|4"T|) on 8 quantities L(r), £(r), t> (r), w(r), fAr(^), 
■tu±(?")> ^( r ) ; p( r ) that are the characteristic turbulent length scale, normalized magnetic 
helicity, matter inflow velocity, turbulent velocity, radial Alfven speed, perpendicular Alfven 
speed, temperature, and density. I rewrite the equations once again in terms of named 
quantities: 

Anpvr 2 = M, (52a) 
vv , , V | R ( P T)> | (p U y r | (ripvl ± )> (rVU _ 
r 2(r — r g ) 2 /i p 3p pr 2 2pr 4 

vv' r + rg ° + w' r + ^uu' r + 2(v 2 A± )' r = with (52c) 
2(r r g ) 3 

RT ( 3ir 3 (0 _1 ) + ^i(6- 1 ) \ 5 2 5RT 5 2 

w — -t(°- 54 9(^(8-1) - 1) +1 - 69 j + r 01 «=^-= 15 -+ g « a , 

t ><r 4 ); = 3.03<.exp(-fl - (0.70< + 0.58(^ ± - ?^>^> _ 064 / RT^a^ 3/2 
pr 4 L Loo V A* / 

(52d) 
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vpr 2 



IlV 3.03v 3 A± exp{-O-(0-^v 2 A± + 0.29{vAr-VA±)vA±)u 0.64 f RT^a^ 3/2 



/WVr ^ ^ \ /J, 



(52e) 

v 2/3 (jf^Y = 1.14^ 3 - 0.09« + 2^>exp(-Q _ L93 / M^ \ 3/2 

+ 2^ ± )£ = 3^00^=^, (52g) 

L^L^M-expT--^^. (52h) 

Here = kT/m e c 2 . Since my prescription for external driving of turbulence is Q + = const, 
I take v p and L to be constant in the source terms. Relativistic wr (eq. |14j ) and non- 
relativistic w^r (eq. [15] ) values of enthalpy w are employed. In the next section I describe 
the values of boundary conditions and parameters for the equations I solve. 



3. BOUNDARY CONDITIONS AND PARAMETERS 

The system (|52|) consists of 5 differential and 3 algebraic equations and should be in- 
tegrated inward from some outer boundary at r x . This requires knowledge of at least eight 
constants. Seven of them are the values "at infinity" L^, T^, p^, C.^, Uoo, tUroo, ^a±oo- The 
eighth is the accretion rate M. It is usually determined by some extra condition and is not 
adjustable. I assume isotropic turbulence with Ek = Em at the outer boundary. Therefore, 

/ -RTooCTocA ^ 2 (ZRT^o^ 1/2 
vatoo = v A ±oo = I I and = I - I , (53) 

and I have one model parameter <7oo instead of 3 velocities tUroo, valoo, and Uoo. Another 
adjustable parameter of the model is 7 that determines the size of energy containing eddies 
L near the object (eq. |52hj ). 

Parameter 7 is not free in principle, but its value cannot be determined within the 
proposed theory. Neither there exist anisotrop ic MHD simulations that could provide 7. All 



simulations to date sh ow 7 to be within 0.2-=-2 (ITennekes fc Lumleylll972uLandau fc Lifshitz 



19871 ; iBiskampI 120031 ) in both HD and MHD case. I assume the same range of 7 in my 



calculations. 



-16- 



3.1. Outer Medium Transition 



Bondi radius 



t'B 



with 



( 



5RT n 



V 3/i 



1/2 



(54) 



is the natural length scale of the spherical accretion flow (IBondil Il952l ). Density p and 
temperature T of plasma are constant for radii r > r^, because gravitation al energy an d 
gas regular kinetic energy are negligible there compared to gas internal energy (jBondilll952l ). 
Averaged magnetic field and averaged random velocity are also constant for r 3> r_g, because 
constant external energy input balances dissipation in this region. As a consequence, £ = £oo 
and L = for r ^> tb- 

I set the outer boundary at r x = 3rg, where matter is almost uniform. Length scale 
should be determined from known external energy input Q + and outer magnetization 

(Too. However, Q + is not known. I assume for simplicity = jtb, so that L changes its 

behavior near tb together with temperature and density. 



Bondi radius is about tb ~ 3 x 10 5 r g for our Galactic Center (IGhez et al.ll2003l ). The 
properties of gas at 3rg are somewhat constrained from observations. I take the values for 
uniformly emitting gas model with temp erature 1.5 x 10 7 K, electron and total number 

densities n eoo = 26cm~ 3 , rioo = 48cm -3 (jBaganoff et al.l 120031 ) at r x = 3rs that corresponds 
to 5" in the sky. The presence of dense cold co mponent can make t he average temperature 



much lower and the average density much higher (ICuadra et ajJl2006J ) , but I am leaving these 
uncertainties for future research. 

Expanding and colliding hyperalfvenic stellar winds provide magnetic field into the 
region. Its strength near Bondi radius is not known. Only the very general estimate can be 
made. Matter magnetization is likely to be lower than the saturation value of — 1. I take 
the values in the range = 0.001 4- 1 to cover all reasonable magnetization states of matter 
at 3r#. If magnetic field is rather a product of decay than dynamo amplification, then the 
local dimensionless helicity £ may be close to unity. I cover the range £oo = 0.001 4 0.5 in 
simulations to determine the possible dynamical significance of non-zero magnetic helicity. 



3.2. Transition to Rotationally Supported Flow 



The system of equati ons (1521) has the same property as spherically symmetric system of 
hydrodynamic equations (jBondilll952l ): subsonic solution exists for all accretion rates M up 
to maximum M*, transonic solution is valid for the only value M*, and no solution exists 
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for M > M*. The solution with 

M = M* (for transonic solution) 



(55) 



is preferable, because it has the highest rate of energy transfer towards the equilibrium state 
of the system matter-SM BH. The same argument is valid for a general hydrodynamic nozzle 
(ILandau &: Lifshitzl 119871 ). It is reasonable to expect that maximum mass flux solution for 
system with magnetic field fl52|) also obeys the condition fl55l) . However, even small amount 
of angular momentum can change the picture. 

Every real astrophysical accretion flow has non-zero specific angular momentum at the 
outer boundary 

/ = Ar g c, or equivalently, / = fi^cir^cir, (56) 

where r cir is a radius where matter becomes rotationally supported and fxdr is Keplerian 
velocity at r c ; r . General Newtonian expression for Keplerian velocity at radius r is 



(57) 



V K = C\ -f- 




At larger radii r > r cir angular momentum exerts relatively small force Fi oc l 2 /r 3 on 
plasma, sinc e Fj decreases with r adius faster than gravitational force F g oc r g c/r 2 . Numerical 
simulations ( ICuadra et al.ll2006l ) suggest r cir ~ 3 x 10 3 r g for our Galactic Center. 



When angular momentum (eq. [55] ) is large, A ^> 1, it should be able to travel outward 
through the outer quasi-spherical solution by means of rip component of stress tensor t a p. 
The angular averaged form of this component is 

< B r B± >q 



An 



(58a) 



where I neglect the kinetic part for the estimate. It can be transformed with the aid of 
Schwartz formula < xy >< V < x 2 >a/< y 2 > into inequality 

B r B i 



tr(f> 



Ait 



(58b) 



with definitions ff2TT) of rms B r and Bj 



Let us take a disk (jShakura fe Sunyaevl Il973l ) with height H and write the angular 
momentum transfer equation as 

d(r 2 Ht r4> ) 



dr 



0. 



The result of integration is (IGammie &: Pophamlll998h 

Ml = AnHr 2 U^ 



(59a) 



(59b) 
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in case of large dimensionless angular momentum A ^> 1 (iGammie fc Pophamlll998l ). I take 
specific angular momentum I from equation fl56l) and the accretion rate to be 



(60) 



M = 2-KrHpv. 



I substitute angular momentum / from relation (1561) . accretion rate M from equation (160|) . 
Alfven speeds from definitions (I24"|) . Keplerian velocity from equation (|57|) . and inequality 
(I58bl) on t r<} (, into angular momentum transfer equation (159bl) to obtain 



vava± 



2X, 



X<1 



(61a) 



that should be valid at any radius r. Sometimes, this inequality is valid for r > r cir if it is 
valid at r C i r , so that condition (16 lap can in some cases be simplified to 



wk 

VaVa± 



< 2 



&t '"cir ■ 



(61b) 



Height of the disk H cancels out of final expression, thus conditions (IBTj) are approximately 
valid even for flows with H pa r. Such flows are likely to describe the realistic transition region 
from outer quasi-spherical inflow to inner rotational solution. There are no extra degrees of 
freedom to put conditions on the surface of compact object, so I consider an object to be 
effectively a black hole. 

Condition of angular momentum transport flBTT) may be stronger than maximum accre- 
tion rate condition fl55|). This depends on the value of specific angular momentum I and 
viscous a parameter (jShakura fc Sunyaevlll973l ). Viscous a is approximately a ~ x ° ac- 
cording to my definitions (1161) and (I61al). If a > 0.5, then acc retion proceeds without direct 
dynamical effect of rotation (INarayan. Kato. fc Honmalll997l ). Thus, two types of solutions 
are possible: 



• maximum accretion rate solutions that describe radial flows with small angular mo- 
mentum I < cr g or large viscosity x cr > 0.5 (subsection 14.11) . 

• flows with the rotational support that work for large angular momentum / 3> cr g and 
small viscosity x ° ^ 0.5 (subsection 14. 2p . 

The condition (1611) gives a crude estimate of the inflow velocity and accretion rate M, since 
it assumes specific angular momentum to be constant down to r cir . As matter travels to r cir , 
the amount of specific angular momentum left becomes smaller. Nevertheless, I calculate 
the solutions with effective angular momentum transport using condition flBTj) to illustrate 
the dependence of accretion rate on model parameters for the rotating flow. 
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4. RESULTS 

4.1. Maximum Rate Solution 

Let me first disregard the angular momentum transport condition (1611) and calculate 
the flow with small angular momentum I <C r g c, when mean rotation is not dynamically 
important. 

The system of equations I solve ( 1521 can be rewritten as 

{FiYr N t (F,r 



D 



for i = 1..8. (62) 



Here Fi(r) are 8 functions I solve for, A^(F, r) are function- and radius- dependent numera- 
tors, and 



v 2 



D = 1 -y-2 ( 63 ) 



is a common denominator. Critical velocity V s is 

5u 2 

V? = c% + 2v 2 A± with c 2 sg = c 2 s + —. (64) 

Effective sound speed c sg is equal to that of plasma with effective particles velocity v 2 g = 
v 2 + u 2 . 

According to the maximum-rate condition ( 1551) I search for a smooth solution that has 
a sonic point at some radius r s . The condition at r s is D(r s ) = 0. Zero denominator requires 
all the numerators A^(F, r) to be zero at r s . It can be shown from system (1521 that all eight 
conditions A^(F(r s ),r s ) = collapse into just one, what indicates that maximum accretion 
rate solution is smooth. Two equalities 

D(r s ) = and iVi(F(r s ), r s ) = (65) 

give the missing 8-th condition on M for system fl52l) and the sonic radius r s . Thus, I have 
7 conditions at the boundary at 3tb and 1 condition somewhere in the region. I employ the 
shooting method to search for M and r s that satisfy the relation (|65l) . 



I obtain the Bondi hydrodynamic model (jBondilll952l ). if I set all Alfven velocities and 



turbulent velocity to zero and use non-relativistic prescription for enthalpy w^r (eq. |52cj ). 
Therefore, the accretion rate M equals Bondi accretion rate of monatomic gas 

M B = \ry PoQ C^-^ ' » 4 x lO^Moyear" 1 (66) 
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in the limiting case of no turbulence. The number is ca lculated for the Black Hole in our 



Galactic Center with ry = 1.1 x 10 12 cm f lGhez et al.ll2003f ). T = 1.5 x 10 7 K, and n « 48cm" 3 



( IBaganoff et al.ll2003l ). Accretion rate M appears to be lower than Mb when turbulent 



energy is non-zero (Fig. [TJ). 

Inhibition of accretion by turbulence has the following explanation. First, energy 
of magnetic field i ncreases inward, therefore it exerts back-reaction force stopping matter 



( 1Shvartsmanlll97ll ). Second, magnetic field serves a very effective mechanism of energy con 



versio n from gravitational to thermal via dissipation of turbulence ( jlgumenshchev fc Narayan 



20021 ). Larger thermal energy corresponds to larger gas pressure that also stops matter. 
Within the deduced model I can estimate the actual decrease of accretion rate M from 
Bondi value Mb- 

I take my reference model to have the values 7 = 1, cr^ = 1, = 0.025 of, correspond- 
ingly, dimensionless scale of turbulence, outer magnetization, and outer magnetic helicity. 
The found accretion rates are 0.14Mb for non-relativistic equation of state and 0.24M# for 
relativistic equation of state. I can now consider the whole ranges of all three parameters 
and explain the observed correlations between them and accretion rate M. 

Larger flow magnetization a results in lower accretion rate M. Larger magnetic field 
and turbulent velocity field exerts larger back-reaction force on matter. Also, transformation 
of gravitational energy into thermal happens more readily if magnetization is larger. Larger 
thermal energy means larger gas pressure and larger back-reaction force on matter striving 
to fall onto the central object. 

Several factors lead to higher magnetization. Larger outer magnetization makes 
magnetization in the entire flow o larger. Then larger dissipation length scale 7 allows 
for smaller dissipation of magnetic field. Larger magnetic helicity £ also lowers magnetic 
energy dissipation and leads to larger magnetization a. These correlations can be observed 
on Figure [TJ Increase of the relative length scale of energy containing eddies 7 from 0.2 to 
2 results (Fig. [Tji) in about 2 times drop in accretion rate M. Accretion rate stays constant 
(Fig. [Tb) at small values of outer magnetic helicity ^ . However, M drops an order of 
magnitude as turbulence approaches highly helical state at outer boundary 3r B with ^ 
close to 0.5. The dependence of M on outer magnetization o"oo is not quite steep: accretion 
rate gradually decreases about 4 times as outer magnetization increases 3 orders of magnitude 
from 0.001 to 1. Surprisingly, accretion rate does not rise to Mb (Fig. [lb) even for very small 
outer magnetization ~ 0.001 for non-relativistic equation of state. Even small outer 
magnetic field increases inwards and influences flow dynamics. 



Accretion rate is systematically about 40% higher (Fig. [TJ for relativistic equation of 
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state (solid line) compared to non-relativistic equation of state (dashed line), because magne- 
tized system has some properties of a non-magnetized one. Formula for Bondi mass accretion 
rate (I6"61) is valid only for non-relativistic monatomic gas that has an adiabatic index T = 5/3 . 
Accretion rate is higher for lower T and is about 3 times larger (IShapiro Teukolskylll983l ) 
in case of ultrarelativistic particles with adiabatic index T = 4/3. Accretion rate M is deter- 
mined by relation (|65|) at a sonic radius r s that is smaller than 10 3 r g (Fig. [TH). Electrons 
become relativistic at somewhat larger radius about 10 3 r g in the solutions of system (152|) . 
This leads to gas adiabatic index T (magnetic field is disregarded) lower than 5/3 at sonic 
point r = r s . Thus accretion rate is considerably larger in case of relativistic equation of 
state. 

It is also instructive to trace the dependence of sonic radius r s on parameters. Sonic 
radius for hydrodynamic accretion of non-relativistic monatomi c gas is equal to several 
Schwarzschild radii r s = 2 4- 10r g (IBeskin fc Pidoprygoral Il995l ). Sonic radius is a con- 
siderable fraction of rs for a gas with a diabatic index V substantially smaller than 5/3 
for non-magnetized accretion (IBondil 119521 ) . Magnetized accretion has the same properties. 
Non-relativistic EOS (solid line) results in very small sonic radius r s = 7-r llr g (Fig. [TH). 
Sonic radius for relativistic EOS (dashed line) is r s = 300 4- 1200r g about the radius where 
electrons become relativistic r ~ 10 3 r g . The value of sonic radius drops several times as 
plasma outer magnetization (Too increases from 0.001 to 1. As outer magnetization in- 
creases, accretion rate drops (Fig. [lb), because density p and gas inflow speed v decrease. 
Then effective sound speed V s equals the inflow speed v at a point closer to the black hole. 

Inflow velocity v as well as other characteristic velocities of the flow are depicted on 



Figure [2] as functions of radius r for the reference model with a c 
All velocities are normalized to the free-fall speed 



1,7 



1, £o 



v ff = c 



r — r„ 



0.025. 



(67) 



I also normalize perpendicular Alfven velocity v a± and turbulent speed u to one dimension. 
Horizontal line on Figure [2] corresponds to radial dependence r -1 / 2 . 

Inflow velocity v monotonically increases inwards, whereas sound speed c s monotonically 
decreases with intersection almost at the sonic point. Radial Alfven velocity VAr, perpen- 
dicular Alfven velocity v a± and turbulent velocity u (Fig. [2]) start out as constants from 
the outer boundary at 3r B , where turbulence is sustained by external pumping. Then these 
velocities increase and deviate from one another. Radial Alfven velocity v Ar appears to be 
much larger th an va i and u in the inner accretion region. This fulfills the expectatio ns of 
earlier models (jShakura &: Sunyaevl Il973l ; IScharlemannl Il983l ; IBeskin fc Karpovl 120051 ). At 
small radius turbulence is driven by freezing-in amplification of magnetic field and random 
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velocity. Left-hand sides of turbulence evolution equations f!52dj) . (I52el) . and f!52fj) dominate 
over corresponding terms with external driving for radius r < 10 4 r g . Internal driving of v& r 
is much more effective than driving of vai_ and u. Therefore radial Alfven velocity va t is 
larger than other two speeds. This refutes any model with isotropic magnetic field. 

Several pairs of lines intersect on velocity plot (Fig. [2]). I consider three main intersection 
points for the reference model with = 1, 7 = 1, ^ = 0.025, and relativistic EOS 
(Fig. [2h). Crossing of inflow velocity v and sound speed c s occurs almost at the sonic point 
at r s , determined by relation 0651) with critical velocity V s (eq. [M])- No plasma waves can 
escape from within the region with high inflow velocity v > V s . Approximately c s ~ V s at 
sonic point r s ^6x 10~ 4 re, because of low magnetization o « 20% in that region (Fig. [3^,). 
Alfven point is determined by equality v = va t at radius r^. Alfven waves cannot escape from 
within the region where inflow speed is greater than radial Alfven speed va t - Equality holds 
at relatively large radius ta ~ 0.03r#. The third combination of the same three velocities 
also gives a characteristic intersection point. Radial Alfven speed v Ar increases faster inwards 
and becomes equal to sound speed c s at about r ps 4r g . Further relative increase of vat leads 
to magnetic energy dominated flow, what can be traced on magnetization plot (Fig. [3h). 

Figure [3^ shows evolution of plasma magnetization o with radius r for the reference 
model. Thermal energy equipartition assumption does not hold, id est turbulent energy 
does not equal to constant fraction of thermal energy o 7^ const. Magnetization a varies 
more than one order in magnitude from 0.07 to 3. It starts out at initial = 1 at 3r^, 
where turbulence is supported by external energy input Q + = const. Then a deviates down 
as r decreases. Magnetization a drops, because length scale L decreases with radius r that 
causes turbulence to decay faster. At about 0.03re magnetization starts to rise as internal 
turbulence driving takes over. Inflow velocity v slightly deviates up from Alfven velocity va 
as r decreases. Since internal driving rate is proportional to v (left-hand sides of equations 
(152dl) . fl52el) . and (152f|) dissipation rate is proportional VA r , parameter a grows slightly with 
decreasing radius. The growth is about a factor of 5 for 3.5 orders of magnitude decrease in 
radius. Magnetization o jumps up in the region very close to the event horizon of the black 
hole. However, this jump may originate from inconsistent treatment of General Relativity. 

The dependence of magnetic helicity £ on radius is shown on Figure [3b. Helicity £ 
behaves almost reciprocally to magnetization a from Figure [3K Such a behavior can be 
seen from magnetic helicity equation (|52gp. Magnetization a decreases order of magnitude 



during the transition from externally supported to internally supported turbulence around 
r « 0.03r#. Magnetic helicity £ also increases an order of magnitude from 0.025 to 0.2. 
Then £ gradually decreases down to initial value. Thus magnetic helicity £ does not change 
dynamics if it is initially small ^ < 0.1. Only when ^ is large, accretion rate drops. 
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Deviation of inflow velocity v from the free-fall scaling r -1 / 2 makes a density profile 
in magnetized flow different from that in standard Advection Dominated Accretion Flow 
(ADAF). I consider the flow where energy is only advected inward. Nevertheless, I obtain 

p oc r~ c with C ~ 1-25 (68) 

almost independently on the parameters or the equation of state, somewhat shallower than 
p oc r" 15 in ADAF. 

The only question left is how well this flow with maximum accretion rate can describe 
the real situation with large angular momentum I. Given the solution of the system fl52l) I 
can check whether the condition for effective angular momentum transport condition (I6T1) 
holds. Condition flBTj) breaks when evaluated for maximum-rate solution with parameters 
£oo ; c"oo ; and 7 within the chosen ranges and circularization radius r cir > r g . This means a 
flow with maximum accretion rate is unable to effectively transport the angular momentum 
outward. The same conclusion can be made simpler. The transport of angular momentum is 
a magnetic process. So, I can be transported only by Alfven waves. However, Alfven waves 
cannot escape from the region within rs 0.03r^ from the compact object that makes 
angular momentum transport impossible even from quite large radius. 



4.2. Solution with Effective Angular Momentum Transport 

Solution with large outer angular momentum I 3> r g c and small viscosity may have prop- 
erties, substantially different from those of maximum-rate solution. The actual details of the 
solution and allowed accretion rate depend on how this angular momentum is transported. 
For the simple estimate I suppose that the accretion rate is determined by the equality in 
angular momentum transport condition flBTj) . Maximum accretion rate M for condition fIBTl) 
appears to be about two orders of magnitude lower than Bondi rate Mb (eq. [66]). 

I add one parameter in modelling: unknown circularization radius r cir for specific angular 
momentum I (eq. [SB]). I take it to be r cir = 10 3 r g for the reference model. Plots of the 
accretion rate verses model parameters are shown on Figure HI Dependencies for the rotating 
solution (Fig.HJ) have the opposite slopes to those for the maximum-rate solution on Figured! 
Accretion rate M increases with increasing outer magnetization o"oo (Fig. Hb) and increasing 
outer magnetic helicity ^ (Fig. Hb). Both effects lead to higher plasma magnetization a. 
I showed in the previous subsection 14.11 that the magnetic field plays an inhibiting role on 
matter inflow, and that the larger the magnetic field is, the smaller the accretion rate M 
is. However, the correlation between the magnetic field and accretion rate is the opposite in 
case of the rotating flow. Accretion rate quantitatively agrees with relation for ADAF flows 
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M ~ aM B ~ crxM B flNaravan. Kato. fc Honmalll997h with a ~ 0.01 at r cir (Fig. W)- 



The allowed by condition (I6T1) inflow speed v is proportional to the product of radial 
Alfven speed va t and perpendicular Alfven speed va±- Larger magnetic field results in larger 
transport of angular momentum outward, so larger inflow velocity v and larger accretion 
rate are possible. Larger outer magnetization o"oo and larger outer magnetic helicity £oo both 
lead to higher magnetization a and higher magnetic field. Inhibiting effect of magnetic field 
is smaller in case of lower accretion rates M and lower inflow velocities v. Lower v results 
in lower relative driving of turbulence that makes magnetic field weaker. Weaker magnetic 
field has weaker influence on dynamics. In sum, larger magnetic field B results in larger 
accretion rate M, when it needs to transfer angular momentum. 

The dependence of M on length scale 7 is obscured by the dependence of external 
driving on 7. Accretion rate M is smaller for smaller magnetic field, but the state of low 
magnetization can be achieved in two different ways. Firstly, magnetic field decays faster 
when L decreases. However, the plasma at circularization radius r c ; r = 10 3 r g is still par- 
tially influenced by the outer boundary conditions. Internal driving does not depend on L, 
whereas external driving is stronger and magnetization a is higher, when L is small. The 
described two effects balance each other and make accretion rate M almost independent of 
dimensionless length scale 7 (Fig. H^). 

Accretion rate M decreases with the decrease of circularization radius r cir (Fig. HU) for 
non-relativistic equation of state. To explain this, I trace on Figure \5jp all the quantities 
that enter angular momentum transport condition (I61b|) for the reference model. Velocities 



normalized by the free-fall speed (eq. |67j ) are shown on Figure [5b. Inflow speed v and radial 
Alfven velocity VAr reach free-fall scaling at about 0.02re. Only perpendicular Alfven velocity 
va± has a different dependence on distance from the central object for r < 0.02r B . Because 
va± decreases with radius, the allowed v and M are smaller for smaller circularization radius. 

However, the accretion rate increases for small circularization radii for 1-T equation of 
state (Fig.HU, solid line). This is the consequence of the decreasing gas adiabatic index, when 
electrons reach relativistic temperature s. Solutions with lower adiabatic index are known to 
have larger accretion rates jBondilll952h that is equivalent to the lower inflow speeds v in the 



solutions for the fixed matter inflow rate. Velocity v (Fig. [5k.) starts deviating down from 
the self-similar r~ 1//2 solution at approximately 10 3 r g , making the solutions with higher M 
possible. In fact, condition (16TT) for the solutions with small r cir becomes critical at some 
fixed point r d > r cir instead of reaching equality at r cir (eq. [5Tb]). Therefore, according to 
condition (loTh ) , maximum value of the inflow speed grows with the decrease of circularization 
radius as v oc r cir , explaining the rise of accretion rate for small r c \ r (Fig. HU, solid line) for 
1-T equation of state. 
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Solution for non-relativistic equation of state, in turn, possess its own feature. Self- 
similar flow (see Appendix [B]) settles in at 10 3 r g , making accretion rate almost independent 
on circularization radius (Fig. SU). Magnetic helicity £ in such a flow is a number about 
unity what is consistent with self-similar solution obtained in Appendix [B] Self-similar flow 
can not establish for 1-T equation of state, because relativistic effects become important 
before it establishes and break self-similarity. 

In fact, magnetization a and magnetic helicity £ (Fig. ED are not constant at small 
radii for correct 1-T EOS, because these relativistic corrections work. At about O.Olre 
magnetization reaches almost constant level a ~ 0.02 (Fig. [6^,) and then starts to slightly 
deviate down, because equilibrium o for matter with lower gas adiabatic index T < 5/3 is 
lower. Magnetic helicity £ behaves (Fig. Eb) the opposite way to magnetization o : magnetic 
helicity reaches £ ~ 1.5 at 0.01r B and starts to slightly deviate up as the radius decreases. 



5. DISCUSSION OF THE MODEL 

I present the sophisticated analytical model to determine the properties of spherical 
magnetized accretion. The common assumptions of magnetic field isotropy and thermal 
equipartition are released, but many assumptions are still left. As usually in fluid dynamics 
a lot of simplifications are made during the course of elaboration. The validity of almost 
everything can be questioned. The system of equations (|52|) may not describe the real flow 
(subsection 15. ip or may have some inaccuracies (subsection 15. 2p . Gas cooling may not be 
neglected (subsections 15. 3p . Convection and diffusion may change the flow structure (subsec- 
tion The equation of state was also found to influence the dynamics (subsection 15. 5p . 
Let me discuss all these topics and determine the practical significance of the model. 



5.1. Real Flow 



Presented model is partially applicable to the real systems. It may describe some gas 
flows onto Supermassive Black Holes in Low Luminosity Galacti c Centers, in particu lar in 
the center of our Galaxy. These flows are geometrically thick (INarayan fc Yil I19951 ) and 



may have low angular momentum (IMoscibrodzka. Das, fc Czernyll2006l ). However, the real 
flows may have properties that my model cannot handle in its current state. First of all, 
the sources of matter and external driving should be explicitly accounted for. Secondly, the 
self-consistent angular momentum transport theory is needed. 



The material is mainly supplied to the central parsec of the Milky Way by stellar winds 
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( Quataertll2004r). The wind-producin g stars have a broken power-law distribution as a func- 
tion o f radius (IBaganoff et al.l 120031 ) . Some stars are as close to the central black hole as 



O.Irs (IGhez et a 



2003 ). The st ars supply too much material to be accreted, therefore there 



exist an outflow (lQuataertll2004j ) . Bondi radius coincides with the radius where inflow starts 
to dominate outflow in numerical simulations with the accretion rate M ~ 10 _6 M©year _1 
(ICuadra et al.ll2006l ). Maximum accretion rate in the solution with zero angular momentum 
is 0.2Mb m lO~ 6 M year _1 and 0.01M# for the rotating flow. So that the transition from 
the outflow to the inflow happens at r > 10 5 r g . 

I can show that outflow from r > 10 5 r g does not change the accretion rate from calcu- 
lated. Outflows substantially alter the value and the sign of inflow velocity v in the system 
(|52p . However, the differences in inflow velocity do not influence any other quantity as long 
as three conditions are satisfied: 



1. v is much smaller than gas particles velocity v p , bulk kinetic energy of gas is negligible 
in the outflow region, 

2. external driving of turbulence Q + dominates over internal driving there, 

3. condition on M is set in the inflow region. 

The first two conditions are satisfied down to r ~ 10 4 r g (Fig. [2] and Fig. [5]). The third 
condition holds for maximum rate solution, because condition on M is set at the sonic point 
about 10 3 r g from the central object. It also hold for the solution with angular momentum 
transport, because the condition on M is usually set at the inner boundary 10 3 -r- 10 4 r g . All 
three above conditions hold, hence outflows of stellar winds do not substantially change the 
accretion rate or any quantity in the system. 



5.2. Treatment of Magnetic Field 



The long history of accretion theory has many accepted models based on ideas, ex- 
tended beyond the area of applicability of these ideas. For example, general relativity 



was substituted with Pa czynski-Wiita gravitational potential (IPaczhynski fc Wiital Il980 



Sha kura fc Sunyaev 



rtNarayan fc Yi 



1995 



973). Magnetic fiel d was long treated similar to the normal matter 



Coker fc Meliall2000l ). Displacement current was n eglected in magneti c 



19731 : lLandau fc Lifshitd 119871 : lLandau. Lifshitz fc Pitaevskiil 11984 : iBiskamd 120031) . Gvroki- 



field dynamics that allowed to treat magnetic field without electric field 
System of viscous equations describe viscosity by a single param e ter ( Sha 



Scharlemann 



19831) 



aira fc Sunyaev 



netics is used to solve the problems with non-Maxwellian distribution functions (jSharma et al 
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2007 ) 
2Q02h 



power-law non-thermal electrons are usually present in plasma (I Yuan. Markoff. &: Falcke 



Described above model is extended in several ways, mainly with regard to magnetic 
field. Isotropic MHD system of turbulent equations (127]) describes the real box collisional 
turbulence quite well, because it corresponds to convergent set of simulations. Collisionality 
assumes that medium behaves like many particles with short-range interactions. However, 
astrophysical medium of interest is always collisionless with prevailing long-range interac- 
tions. I inconsistently use the results of numerical simulations of collisional MHD (eqs. [HE]) 
with magnetic resistivity vu on the order of viscosity u, because the realistic simulations of 
collisionless plasma turbu lence are not done and are unlikely to be done in the near future 
rtSchekochihin et all 12004 ). 



Observations of astrophysical turbulence may give more information than numerical 
simulations. A special case of collisionless plasma is plasma with random kinetic energy 
much smaller than random mag netic energy. This regime is a good picture of Sun corona 
with all plasma effects into play (1 Aschwandenl 120051 ) . Dissipation of magnetic loops with low 
kinetic energy proceeds mainly via reconnections. The timescale of reconnective dissipation 
was found to be 

r rec « 20— (69) 
va 



in solar flares (Noglik. Walsh, fc Ireland! 120051 ). The same number was also predicted by 
Lazarian fc Visnniad ( 19991 ). Collisional MHD turbulence has much smaller dissipation 
timescale 

Tdiss ~ 1— (70) 

va 

(eqs. [271 EH]). Plasma has large kinetic energy in the outer region of accretion flow, where 
turbulence is externally supported. Timescale Tdi SS (eq. [70] ) may be appropriate there. 
Kinetic energy Er decreases to smaller radii and magnetization o increases (Fig. [2]) in 
case of zero angula r momentum (14.11) . Accretion flow there may resembles solar Corona 
( lAschwandenl 120051 ). Dissipation timescale may increase order of magnitude and be close 
to r rec (eq. [ni])- This increase would lead to much lower accretion rate, because higher 
magnetic field leads to lower M. Matter infall may eventually proceed through channels of 
lower magnetic field (jlgumenshchevl 120061 ) . 



Even if I assume that box isotropic turbulent system of equations ( 1271) with coefficients 
( 155]) is applicable to isotropic turbulence, there are at least four complications in building 
the full anisotropic theory. 

First of all, I need to introduce arbitrary coefficients c U B2, cbb2, cb U 2 to describe 
isotropization of anisotropic magnetic field and anisotropic energy transfer between mag- 
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netic field and fluid motions. Reasonable values of these coefficients were taken to satisfy 
rather loose analytical tests (Appendix IAJ . However, changes in these coefficients do not 
lead to dramatically different accretion rate or flow structure. Setting cbb2 = cbbi instead 
of cbb2 — leads to only 10% of M change for the reference model. All seven introduced 
coefficients c xx may themselves depe nd on anisotropy of the magnetic field. T he details of 
anisotropic MHD are still debatable (IGoldreich fc Sridharlll995l ; lBoldyrev II2006I ). I leave the 
incorporation of anisotropic MHD model into accretion theory for future work. 

Secondly, the presented theory is not general relativistic. Accretion rate M appears to 
be insensitive to the choice of gravitational potential. The condition on M is set at about 
10 3 r g in case of relativistic EOS and zero angular momentum I. Sonic point is situated 
close to the black hole at r s = 5 -f- 10r g for non-relativistic equation of state. But 1% 
increase of M leads to the sonic point at r s > 100r g , independent of the way to mimic 
general relativity. However, the region near the black hole is important, because part of 
synch rotron IR radiation as well as part of radio emission comes f rom several Schwarzschild 
radii (INarayan et al.lll998l ; iFalcke fc Markoffll2000l ; iMarrond 120071 ). Thus, to fully constrain 
theory by observations general relativistic magnetohydrodynamics is a must. 

In third, magnetic helicity H involves numerous complications. Magnetic helicity evolves 
in the region that is frozen into matter. The distance Lu between radial boundaries of this 
region is proportional to inflow velocity v, thus Lu increases with increasing v and at some 
point Li I > r, whereas size in the angular direction is about L = 77*. A part of the region 
is getting sucked into the black hole, while a part is still situated at fairly large radius r. 
Equation of magnetic helicity evolution ( 52g ) holds only if I assume even redistribution of 
magnetic helicity over the mass of plasma. This holds for frozen magnetic field, but in reality 
diffusion and convection are present. Diffusion may change the results for H (eq. |52g| ) as 
well as for the entire flow pattern. I also leave these uncertainties for future research. 



In fourth, it was recently suggested by iBeskin &: Karpovl (120051 ) that ions and electrons 

'his is the opposit e 



should be viewed in accretion as confined by magnetic field lines, 
of standard picture where magnetic field lines are frozen into matter (IScharleman 



The former case has higher heating rate of matter under contraction (IBeskin fc Karpov 



o pposit e 
rJll983h . 



2005h. because of conser vation of the first adiabatic invariant I = 3cp^/(2eB) = const 



( jLandau fc Lifshitzlll975l ). Here pt is a particles momentum in the direction perpendicular 
to B. However, only highly magnetized flows with magnetization a > 1 conserve I. Non-linear 
collective interactions of particles in low-u plasma are likely to isotropize their distribution. 
When particles are heated isotropically under contractio n, general Magneto-Hydrodynamics 
(eqs. [HE]) works (jLandau. Lifshitz fc Pitaevskiil Il984j ) and heating rate stays unchanged. 
Magnetization in computed models is below unity (Fig. [3^ and Fig. [6^). Thus application 
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of first adiabatic invariant conservation to magnetized accretion flow seems irrelevant. 

Finally, mean rotation of the flow also creates anisotropy. Because the inner gas rotates 
faster than the outer, MagnetoRotational Instability (MRI) works. It prod uces the additional 



drivin g of magnetic field that may be concurrent to other sources. MRI ( jHawley fc Balbus 



20021 ) has a timescale 



x M , = -i3^) . (71) 



dr 

When MRI timescale becomes larger then dynamic timescale T^ yn = r/v, field amplification 
occurs mainly because of regular shear tangential motion, instead of regular radial motion. 
MRI may be crucial even in the region without rotational support. Full consideration of 
effects of angular momentum on the flow is the subject of the next study. 



5.3. Radiative Cooling 



The system of equations (152]) describes the accretion flow, where all the energy is stored 
in the same piece of matter where it initially was. There is no energy loss by diffusive or 
radiative cooling. But whether such a model is realistic. 

Let me estimate the radiative cooling first. Line cooling is more effective than bremsstrahlung 
cooling for te mperatures about 1.5 x lO 7 ^. Line cooling function is A « 6x 10 _23 n 2 (T/10 7 i^) -0 ' 7 

erg cm -3 s -1 ([Sutherland &: Dopitalll993l ). Thus characteristic cooling time r coo i is 



^cool 



3RTp 
2Afi 



1 x 10 12 s 



(72) 



for our Galactic Center accretion. The dynamic timescale r^ yn = r/v for accretion with rate 
M = 0.1M B (eq. [66]) is 



T d yn = ^~ ~ 5 x 10 10 s 
M 



(73) 

with continuity equation (jUJ) at radius r — tb (eq. [54]). Cooling time is about 20 times 
larger than inflow time in the region where outflows dominate. Nevert heless, anisotropy o f 
stellar winds may lead to significant cooling of some clumps of matter (jCuadra et al. 2005 ). 
Even the disk may form (jCuadra et al.ll2006l ). Careful calculation with line cooling is yet to 
be done. 
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5.4. Convection &c Diffusion 



The system (152]) does not include diffusive or convective transport of quantities. Thus 
the system represents Advection-Dominated flow, where magnetic field and gas can exchange 
energy between each other. The exact model would include transport of momentum, energy, 
magnetic field, magnetic helicity that may or may not influence the dynamics. 

First or all, any type of convective or diffusive motion would happen at a speed v c not 
exceeding the maximum of turbulent speeds, radial Alfven speed v c < vat- This leads to the 
transition from convection dominated to advection dominated flow at several dozens r g in 



the case wit h rotation (lAbramowicz et al.ll2002l ). Correspondingly, inflow speed v becomes 
large v c ~ v ( Gammie fe Popham 19981 ). Transport becomes ineffective at r < r^, where ta 
is the radius of Alfven point. According to Fig. [2^, Alfven point in my spherical solutions lies 
at ta ~ 0.03r#. Thus diffusion and convection are st rongly suppressed in th e inner flow. By 
the same reason, magneto-thermal instability (MTI) (IParrish fe Stondl2005l ) is not supposed 
to play any role for spherical inflow, but may play a role in a case with rotation. For the 
non-conductive convective stability criterion see Appendix ICl 

However, speed of electrons v e may overcome the spee d of sound c, s , so electron con- 
duction may in principle transport energy from within ta (j Johnson fe Quataertl 120071 ). It 
is yet unclear whether electron conduction is suppressed at high inflow velocity v > vat, 
because electrons may be bound to the field lines of tangled magnetic field. The efficiency 
of conduction is a free parameter. If efficiency is close to maximum and conduction is not 
i nhibited, then accretion ra te may be 1 -j- 2 orders of magnitude lower than Bondi rate Mb 
( I Johnson fe Quataertll2007l ). thus accretion rate would be limited by conduction and not by 
backreaction of the magnetic field. Other types of energy transport (IParrish fe Stond 120051 ) 
may kick in for lower accretion rates. The correct calculation with magnetic field and better 
prescription for conductivity is yet to be done. 



5.5. Equation of State 

The difference in accretion rate M between one-temperature relativistic and 1-T non- 
relativistic EOSs is up to 40% for maximum-rate solution (subsection 14.11) and up to several 
times for solution with effective angular momentum transpor t (subsection I4.2p. Solutio n 



with smaller gas adiabatic index T has larger accretion rate M ( jShapiro fe Teukolskylll983l ). 
Gas adiabatic index gradually falls from T = 5/3 to T = 1.43 in case of relativistic EOS as 
matter approaches the black hole. 



However, the electron temperature T e is unlikely to be equal to ion temperature Tj. 
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Electron temperature T e is usually modelled to be lower than (INarayan &: Yilll995l ). This 
two-temperature model has lower gas pressure support and larger gas adiabatic index T than 
1-T model with T = Tj. Lower gas pressure leads to higher accretion rate, larger T leads 
to lower accretion rate. The combination of these two effects is expected to change the 
accretion rate by about the same 40% as between relativistic and non-relativistic 1-T EOSs. 
The exact details depend on the two-temperature model chosen. 



6. OBSERVATIONS 



Proposed quasi-spherical magnetized accretion model is aimed to explain plasma flow 
onto SuperMassive Black Hole Sgr A* in our Galactic Center. Many observations of this 
source are made. These observations reasonably agree with the results of my model. 

A common misconception about Chandra X-Ray observations of Sgr A* exists in liter- 
ature. X-Rays mainly originate in the region that lies further than Bondi radius r# from 
the central object. Thus characteristic density and temperature far from the Black 
Hole can be found (IBaganoff et al.ll2003l ). If one knows the mass M, this automatically gives 
Bondi accretion rate Mb (eq. [66] ) . However, a ccretion rate is no t necessarily determined 
by this formula (IST)]) . unlike some papers suggest (IBower et al.ll2005l ). In my model accretion 
rate M is independent on radius and is smaller than Mb- 



IR (lEckart et al.l 120061 ) and Radio (jShenl 120061 ) observations are difficult to interpret, 
because fluxes in these diapasons depend strongly on the accretion model. Density of mat- 
ter p is better constrained by observations than accretion rate M. The general agreement 
( jYuan. Markoff. fc Falckei 120021 ) is that density p should be low er than in Bondi solution p B 
in th e region close to the black hole. Sol utions with outflows (|Yuan. Quataert. fc Narayan 
20031 ) and Convectively-Dominated flows (IQuataert fc Gruzinovl I200CH ) were invented to ex- 
plain this lower density. Magnetized solution without angular momentum does well the same 
job. Let me consider the reference magnetized model with = 1, 7 = 1, = 0.025, / = 0, 
1-T relativistic equation of state. The ratio of density in a reference magnetized model to 
density in a non-magnetized solution is 



-"magn 



0.27 at 10r P 



P 



(74) 



nonmagn 



Density in a magnetized model is much lower than in a non-magneti zed one. However, all 
types of models can be made to fit the data by adjusting temperature (IQuataert &: Gruzinov 
20001 ). whether advection dominated or convection or outflow dominated. 



Faraday rotation of submillimeter radiation offers a good differentiation mechanism be- 
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tween ADAF flows and flows with outflows or convection. Rotation measure is pro portional 
to bo th magnetic field and electron density and has a relativistic temperature factor (IMarrone 
20071 ). Model B predict s magnetization a = 0.7 and number density n = 2 • 10 7 cm -3 at 3r g 
that is consistent with ( Hawlev Sz Balbusl 2002 ). The observed Faraday rotation measure 
is RM = —6 • 10 _5 rad m" _1 . JMarrone![2007] l. Fitting the relati vistic rotation measu re for 
temperature gives T e = 4 • 10 10 K in excellent agreement with (ISharma et al.ll2007l ). Ac- 
cre tion rate in the refe rence model is about 9 • lO _7 M year _1 , what is 30 times lower than 



in (ISharma et a 



20071) . However, the electron density in my model is close to that in the 



rotating model (ISharma et al.ll2007l ). because inflow velocity in the rotating model is a times 
lower. For densities to agr ee I need a ~ 0.03 that is somewhat smaller than found in numer- 
ical simulations a > 0.2 (IHawley &: Balbusl 120021 ) . This means my solution overestimates 



densi ty n by about a factor of 5, what results in larger t hen observed IR flux (jEckart et al. 
2006 ). Effec ts of angular momentum tr ansport, outflows ( Yuan. Markoff. &: Falcke 20021 ) 



or 



conduction (IJohnson fc Quataertl 120071 ) must come into play to allow for successful fitting 
for both IR flux and Faraday rotation measure. 



7. CONCLUSIONS 

Though many ways of dealing with inefficient accretion were invented, my approach is 
substantially different from all previous efforts. I elaborated the model that 

• has very few free parameters, 

• self-consistently includes averaged turbulence, combining geometrical effects of freezing- 
in amplification with dissipation, 

• ties evolution of random magnetic field and random velocity field to numerical simu- 
lations, 

• connects outer externally supported turbulence to inner self-sustained turbulence, 

• predicts the accretion rates M and flow patterns for the flows with negligible angular 
momentum, 

• gives the order of magnitude estimate of M for large angular momentum flows. 
The model predicts 

• accretion rate M of magnetized fluid 0.2 0.7 of Bondi rate Mb even for small outer 
magnetization a^, 
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• subequipartition magnetic field in the outer part of the flow and superequipartition in 
the inner part, 

• several times lower density than in Bondi model near the central object, what with 
addition of other effects would explain the observations of Sgr A*, 

• half an order of magnitude effect of different equations of state on the accretion rate, 

• unimportance of magnetic helicity conservation, 

• ineffectiveness of convection. Convection and diffusion should be accounted for to- 
gether. 

The next version of the model will include 

• more anisotropic effects, in particular, magneto-rotational instability, 

• two-temperature equations of state, 

• full treatment of angular momentum transport, 

• diffusion of momentum, heat and magnetic field. 
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Let me consider my model in anisotropic incompressible case of box turbulence. I 
substitute —vd/dr = d/dt in equations (122"]) . (123]) . (126]) and set r = const. The box has 
infinite volume. I express some of unknown c xx in terms of known c xx from equations (128]) . 
The system now reads 



A. ANALYTICAL TESTS 




(Alb) 



(Ala) 



This preprint was prepared with the AAS IATgX macros v5.2. 
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d(u 2 ) _ (c uB (v 2 Ar + 2v A±) - C u B2(vAr ~ V Al) 2 )u - C UU U 3 

IT ~ l • (Alc) 

I need to determine three coefficients cbb2, c U B2, and cb U 2 and prove the entire system (I All) 
makes sense. 

There are three kinds of analytical tests divided by the degree of their certainty. The 
tests from the first group have solid physical grounds. The tests from the second group 
represent how turbulence is believed to work, these are the general relations with clear 
physical insight. The third group of tests consists of the order of magnitude relations and 
the disputable ideas. 

The tests of the first group are proven to work. Only one test of this kind can be 
applied to our system. This is the energy decay test. Free incompressible MHD turbulence 
has decreasing with time total energy, beca use energy decrease corresponds to t he increase 
of entropy of the system gas/magnetic field (ILandau. Lifshitz fc Pitaevskiil Il984j ) . 



d_ ( v\ T + 2v\ , + u 2 
dt 



< for at least one of VA r , v A ±, u non — zero. (A2) 



I take sum with proper coefficients of the right-hand sides of system flAlj) . Then I maximize 
it with respect to va_l/vat and Va/u. I find that when 

2c B u2 + c uB 2 > -2.2, (A3) 

total energy decreases with time for any non-zero v An va±, and u. Let me remind the reader 
that all these velocity are non-negative according to definitions ff2^1) . Condition (]A3j) is weak. 
Some tests from the second and the third categories constrain c U B2 and cb U 2 better, thus 
making equation (IA3I) valid. 

The typical test of the second category deals with dynamo amplification of anisotropic 
field. Dynamo action not only amplifies magnetic field, but also isotropizes it. I take 
isotropization condition to be 

^ " VA± \ < 0. (A4) 
dt(v Ar - v A ±) 

Taking expressions for derivatives from system (1A1I) I arrive at 

(c Bu ~ 3c Bu2 )w - VSc B B2(vAr + v A ±) + C BB2 (2VAr + v A ±) ^ (A5a) 

This condition should hold when any speed in inequality ( lA5aj) is much larger then two 
others. Therefore, inequality (1A5aj) is equivalent to 

a/3 

CBu < 3C BU 2, C B B2 < —^-cbb- (A5b) 
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Another second category dynamo test states that magnetic field should always increase, 
if dynamo operates without dissipation or any energy transfer. This occurs when Alfven 
speeds are much smaller than turbulent velocity field u. Positive amplification condition 
then reads 

at v Ar at v A± 

Taking the expressions for derivatives from system (lAll) and applying the limit v A r "C u and 
va± C«I obtain that inequalities (1A6I) are valid for any balance between vat and va± when 

c Bu > 2c Bu 2- (A7) 

Inequalities (lA5b|) and (1A7j) give tight constrains on cg U 2. 

The similar test exists for the random velocity Magnetic field is supposed to increase 
the turbulent velocity in the limit Va t ~ va± ^ u - The correspondent condition 

d fu 2 \ 

— ( — J > for u Ar ~ va± > u (A8) 

reduces for system (1A10 to the condition of constant positive acceleration that initially steady 
magnetic field applies to matter. Finally 

C U B2 < c uB - (A9) 



Decay of isotropic MHD turbulence offers the following test of the second kind. Nu- 
merical simulations show equality of magnetic field dissipation rate and random velocity 
dissipation rate (133]) when initial magnetic energy equals initial kinetic energy. However, 
this equality should be stable, otherwise kinetic and magnetic energy would diverge from 
each other after any perturbation and equality of u and v a would not have been observed. 
Stability condition is 

d(v 2 Ar + 2v 2 AL - u 2 " 



dt (v 2 Ar + 2v AL - u 2 ) 

for VAr = VAX = u. 



< (A10) 



The are no more proven or justified assumptions I can make. I need to make use of 
inequalities (1A3I) . (1A5bj) . (1A7I) . (1A9j) . and (1A10I) and apply unjustified tests. I take the value 



of cbu2 to be in the middle of the allowed interval 



CBu2 = I ( I + \ ) C Bu « 0.29. (All) 



2 \2 3 
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The value of c u b is small compared to the values of other coefficients. There is no physical 
sense in the sharp increase of u 2 build-up when magnetic field becomes anisotropic that would 
be the case for c U B2 <S (—c u b) Turbulent velocity may be expected to increase regardless of 
the direction of magnetic field in equation (j Aloft . This idea leads to |c U £ 2 | < c u b- I take 

C U B2 = (A12) 
for the simple estimate. Similar estimate allows me to set 

CBB2 = 0. (A13) 
In this case isotropization of magnetic field has a timescale about the dissipation timescale. 



B. SELF-SIMILAR SOLUTION 

Let me describe the self-similar solution, when the differential system of equations 
f l52|) can be reduced to the algebraic system. I set the proper scalings of quantities with 
radius and make weak additional assumptions. I introduce the standard dimensionless 
variables T(x), p(x), L(x), aa(x), bb(x),pp(x), vel(x) to replace, respectively, T(r), p(r), L(r), 
u(r),VA, r (r),VA±(r),v(r) as follows: 

(2RT{x)\ 1/2 

T(r) = TooTix), v(r) = velix) I ) , L(r) = (r/x)L(x), 

\ V J 

(Bl) 

VA±{r) = pp{x) I I 

Radius is normalized to Bondi radius (eq. [54] ) as r = tb x. The natural power-law radial 
dependencies of these quantities ( IB II) 

T(x) = T ss x-\ vel(x) = v ss x~ 1/2 , L{x) = 7 , (B2) 

aa(x) = u S sx~ 1/2 , bb(x) = v Ar ssx~ 1/2 , pp{x) = v A1 _ ss x~ 1/2 

make my system of equations (|52[) independent of x under the following restrictions: 

• gravity is Newtonian, 

• external turbulence driving is negligible, 



u(r) 



aa[x) 



2RT(x) 



v Arij) = bb(x) 



( 2RT(x) 
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• equation of state is non-relativistic. 



These assumptions are valid in the intermediate region 10 r g < r < 0.1r B . Gravity is 
Newtonian for r 3> r g . Turbulence driving is mainly internal for r < O.lre (see subsections 
(14. ip . (14.21) and Fig. [Tb, Fig. Eb). Electrons become relativistic at around 10 3 r g . The found 
range of r where all above assumptions hold is small. I can instead consider a non-relativistic 
equation of state with w = wnr (eq. |52cj ) everywhere. This makes standard self-similar 
solution possible from O.lrg down to several Schwarzschild radii r g . 

Dimensionless magnetic helicity £ appears to be constant in self-similar regime. Rela- 
tions (l52hj) . fl521, and flU} lead to 



^ ~ AT (n? _j_ On, 2 ( B3 ) 
^SSKVArSs + ^A±SS) 

Continuity equation (I52ap can be used to obtain the scaling of density p ~ x~ 3 / 2 . Heat 
balance equation (fTT|) reduces to the equality of radial and total perpendicular magnetic 
fields 

viss = 2v 2 A±ss . (B4) 
Euler equation (152b|) gives the formula for self-similar temperature 



T ss = 5/(15 + 10u 2 ss + 9v 2 ArSS + 6v 2 A±ss + 6*4). (B5) 

Turbulence evolution equations (122"j) . (1251) . (126"j) are now treated without source terms. They 
give, correspondingly, three relations 

ZussVAr-ssCBuii - 2v 2 ArSS c B Bi\ exp(-£) + 4u S sCbu22Va±ss + 3v ArS sVssl = 0, 
2u ss (vArSsCBu22 + (c B «ii + cb u22 )va±ss) - v a ±ss{2c B bu exp(-£)v A±5 . 5 + 3w 5S 7) = 0, (B6) 

-m| s c uu + c uB11 exp(-0«-ss + 2 ^a±ss) = °> 

where definitions of Alfven and turbulent velocities (1241) are used. 



Let me first set magnetic helicity to zero £ = and consider four equations (IB4I) and 
( 1B6I) on four velocities wss, v Ar $s, The only solution of this system has all the 

velocities identical zeroes. No self-similar solution is possible for zero magnetic helicity £. 

However, the non-linear algebraic system of equations on £ and velocities (IB3[) . (IB4j) . 
( IB6I) possesses a non-trivial self-similar solution. For the full system ( 15 2 p I need the additional 
condition to determine the accretion rate and solve for radial dependencies of quantities. 
This condition is either condition for maximum accretion rate ( 1641) or condition for effective 
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angular momentum transport ([61]) . I can transform both into self-similar form. Maximum 
M condition fl64|) reads 

5 + 104 5 + 12v a _l S5 = 6^| 5 . (B7a) 
Effective angular momentum transport condition (IBTj) gives 

^ ^ < 1 (B7b) 



4vatSS va±ssVTss 
regardless of circularization radius r cir . 

Let me first find the self-similar solution in case of large angular momentum. I solve 
equality in relation (lB7bl) and 5 equations flB3j) . flB4j) . (1B5I) . (IB6I) for 7 quantities £, T 5S , 



M 5S) v ArSS, v a±ssi 7 W S5 ; an d the product cr^oo. I normalize the results to free-fall velocity 
(eq. |67j ) to be able to directly compare with the numbers on Figure [5b: 

^ = 0.58, ^ (r) = 0.0094, ^ = 0.041, (B8) 
Vf f {r) V3v//(r) *W) 

= 0-029, = 0.0033, (Too^oo = 0.00718 (B9) 

w //( r ) v //( r ) 

for r ^> r g . Figure [5b shows profiles of velocities for the reference model with = 1, 
£oo = 0.025, 7 = 1. The actual velocities on the inner boundary at r = 3 x 10~ 4 rs = 90r g 
are 

,Jr) o.58, M(r) = 0.0033, ^fl= 0.076, (B10) 



//('•) ""' V&v f f(r) "' " '' ;/('•) 

; ( = 0.024, — Vv = 0.0051. 

v ff (r) v ff (r) 

The reference model has cXoo^oo = 0.025 about 3 times larger than in self-similar solution 
(IB8I) . magnetic field in the reference model is stronger. Therefore, higher values of all 
characteristic velocities are expected in the actual solution fIBlOl) . I obtain inflow velocity 



v and radial Alfven speed VAr correspondingly 1.5 and 1.8 times higher for solution (IBlOj) . 
Sonic speeds are the same in self-similar (1B8|) and actual (IB10I) solutions, because almost all 
gravitational energy goes into thermal energy in both cases. However, perpendicular Alfven 
velocity va± and turbulent velocity u do not qualitatively agree with self-similar solution. 
They are correspondingly 1.2 and 2.8 times lower in the actual solution (1B10I) . The naive 
estimate for accretion rate is 

47rp 00 u(r B )ri « 0.05Mb- (Bll) 

This appears to be 8 times larger than the actual accretion rate 0. 0061Mb. Velocity near 
Bondi radius (eq. [54]) is much smaller than self-similar value, what leads to an overestimate 
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of M. Thus, self-similar solution can give an order of magnitude estimates for all character- 
istic velocities of the flow and even for accretion rate M. However, self-similar solution has 
only 2 free parameters instead of 3, because cr^oo is treated as one constant. Therefore, 
solution of the full system (J52J) is required to probe the entire parameter space and to achieve 
more precise results. 

Self-similar solution in case of maximum rate flow with condition (lB7ap does not exist. 
The formal solution of equations (IB7ap . (lB3p . (IB40 . (IB6[) leads to negative product (TooCoo- 
The absence of self-similar solution in this case is reasonable, since the actual solution does 
not exhibit self-similar scalings (Fig. [Tb). 



C. CONVECTION 

Let me elaborate the stability criterion against convection in my model. As I noted 



in the main text (subsection 15.41) . small scale p erturbations o 



by diffusion. Thus high-frequency analysis by IScharlemannl (119831 ) is not appropriate to 



quan tities are smeared out 



determine the convective stability. Timescale of diffusion t^is is 

Tdiff ~ -, (CI) 

u 

where / is the scale of perturbation. As h decreases, diffusion time also decreases and becomes 
smaller than perturbation growth timescale r grow . If Tdiff < r g row, convection is ineffective that 
is likely to happen at small scales h. Thus I need to consider the motion of the large blobs 
of the size h ~ L. 

I consider a blob of plasma displaced at some small Ar from its equilibrium position 
(Fig. [7j). The density of the blob itself changes by Ap blob , when it is moved. The density of 
outer medium changes by Apfl ui d between two positions of the blob. The goal is to calculate 
the difference in density differences Apfl ui d — Ap b i ob between the outer medium and the blob. 
Positive difference Ap fluid — Ap b i ob > for positive Ar > implies convective instability. 
Rising blob of gas is rarified compared to the fluid and buoyant. The results for Ap may be 
affected by external driving that is somewhat artificial in my model. Thus I need to calculate 
Ap in the inner accretion region where external driving is not important. Motion of the blob 
is adiabatic and governed by the same adiabatic dynamical equations (152bl) and (I52c| . as 
the rest of the fluid. I neglect energy, associated with gas regular velocity v. Term v 2 cannot 
be neglected onl y in the region, wher e v approaches sound speed c s . However, convection 



ceases if v ~ c s (INarayan et. al.l 120021 ). I denote by index A physical quantities in the blob 



and by index F quantities in the rest of the fluid. 
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Euler equation f!52bj) results in the following equations on differences in the blob 

^A A (pT) + l -A A (pu 2 ) + ±A A (r 2 pv 2 AA _) - ^4 A A (rV£ ) = (G2a) 
and in the fluid 

jA F (pT) + \a f ( P u 2 ) + ±A F (r 2 pv 2 A± ) - ±A F (rW A ) = 0- (G2b) 

In both equations I take variations between quantities at r + Ar and r. I introduce the 
difference operator 

A() = A F () - A A Q (C3) 

and calculate the variations of all quantities between the fluid and the blob. Subtracting 
equation flC2b ) from equation flC2a ). I find the radial pressure balance in the first order in 
Ar 

-A(pT) + l -A{pu 2 ) + A(pv 2 A± ) - l -A{pv 2 A ) = 0. (C4) 
p 6 I 

Blob of plasma should be in equilibrium also in perpendicular direction, not only in 
radial direction. I use the same technique to deduce it, as I used to derive the radial force 
equation (fi"3l) from general momentum equation ([3]). Component 9 of magnetic force in 
equations ([2]) and ([3]) reads F e = [B x [V x B]] e /(4vrp). I subtract B g (V x B)/(4vrp) from 
it and average over <f> direction. I obtain 

F e = ^SA (C5) 
an pr 

for Bq = B\ and B r Bg = on average over 0. The final form of force balance in 9 direction 
is 

d (R „ 1 „ 1 



Perpendicular force balance flC6j) has the same form in any direction perpendicular to the 
radial vector owing to the symmetry of the problem. I apply operator A (eq. |C3| ) to the 
integral form of perpendicular pressure balance and get 

^A(pT) + l -A(pu 2 ) + ~A(pv 2 A ) = 0. (C7) 

Heat balance equation (fTTI) gives the third relation 



„ V 2 AT - f T ) + V" A " - 17 ) + " A I p ) + = °- < a 
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Expansion or contraction of a blob is non-uniform. Perpendicular b and parallel a sizes 
(Fig. [7]) deform in different ways. Continuity equation for the fluid ([9]) can be written as 

A F p Apv Ar . 
— — + — — + 2 — = C9a 

p v r 

I consider the parcel with constant mass m = pV. Therefore 

^ + ^ + 2 ^ = (C9b) 
p a b 

is the continuity relation for the parcel. Finally I subtract equation (lC9ap from equation 
IC9bl and obtain 

Ap Av Ar Aa Ab . 
— + — + 2 2— =0 CIO 

p v r a b 



for the change of density according to definition (1C3I) . Inflow velocity v is clearly associated 
with the fluid, but I omit subscript F at v . I also omit subscript A at dimensions of the 
blob. 

Now I need to quantify the variation of the turbulent magnetic field and the random 
velocity. I assume that the blob moves at a speed V(r) much higher than the inflow ve- 
locity V(r) ^> v(r), therefore magnetic field does not dissipate in the parcel. Differences of 
turbulence evolution equations fl52d.fl . (I52el) . and (152f[) are 

2 Ap Ar 

2uAu - -u 2 — = — (c uu u 3 - c uBU (v 2 A + 2vA ± )uexp(-g)) (Clla) 
6 p vL 

2^ , A 2 ( Ar Ab\ _ pAr 3 2 



MP v a)+^P v A r = r~( C BBllV A exp(-£)-(c B ullV Ar + 2CBu22VArVAx)u) (Cllb) 

\ r b J vL 

A(pv A± ) + 2pv A± + % - % - f) = (Cllc) 

= ^ 1 ( C BBH t; i± ex P( _ ~ ((cBull + C B u22)v A± - C Bu22 VArVA±)u). 



Magnetic helicity variation does not directly influence the dynamics of the blob. Solving 
the system of 7 equations (jC3j) . flCTD . (jCSj) . flCTUj) . flCTTk bc) on 7 quantities AT, Ap, Av A , Av A± , 
Au, Aa, Ab, I obtain 

Ap correc t _ 2.02 exp(-£)v Ar v A± (v Ar + 2v A1 _) - u(0.39(v Ar v A ± + v\j_) + ^r(l-21^i_L - 0.63m 2 )) 
pAr ~ Ar c 2 s Lv(v 2 Ar + v 2 A± ) 

(C12) 

The actual expression is much longer. I take only the largest terms in the numerator and 
the denominator. 
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Let me compare this result (eq. |C12j ) w ith the naive estimate, when magneti c field 
dissipation increases gas internal energy only (jBisnovatyi-Kogan &: Ruzmaikinl Il974l ). and 
gas pressure balance is used instead of parallel and perpendicular pressure balances (|C4p . 
0C7p . Gas pressure balance is 

A(pT) = 0. (C13) 
Naive heat balance ffl6l) for the unit mass is 

— ( ^-pAT-TAp] « ^(0Alv Ar u+lAQv Ar uv A± +lAuv A± -3m(v Ar +2v A± )exp(-£)-l.Uu 3 ) 

(C14) 

Eliminating AT from relations ( 1C13I) and ( 1C14I) . I find 

1 Apnaive _ 0-61 (v\ r + 2v\ ± ) exp(-Q + 0.23m 3 - Q.%2v 2 Ar u - Q.23v Ar uv Ai _ - 0.28w^ ± 
p Ar c 2 s Lv 

(C15) 

I evaluate the convective derivatives of density (1C12j) and (1C15|) in the inner region of 
the reference solution with angular momentum transport (subsection 14.21) . Parameters of 
the reference model are ^ = 0.025, = 1, 7 = 1, non-relativistic EOS. Correspondent 
velocities are shown on Figure [5b. I take the values (IBlOp of velocities and magnetic helicity 
on the inner boundary of integration at r = 3 x 10 _4 r^ « 90r g . Change of density appears to 
be negative Ap < for Ar > in the result of full calculation (eq. [C12] ). Naive calculation 
shows positive Ap > for Ar > 0. 

^ Pcorrect « -0.2. (C16) 

Naive calculation suggests that the flow is convectively unstable, whereas the full calculation 
under reasonable assumptions indicates a convectively stable flow. 

The calculated result (lC16j) is applicable only to the inner regions of solution with 
angular momentum transport (subsection |4.2p . Excluded external driving is important in 
the outer regions. In turn, solution with maximum accretion rate has large inflow velocity 
v that approaches gas sound speed c s , and convection is suppressed (subsection 15. 4p . As 
a bottom line, either flow appears to be convectively stable on average or convection is 
suppressed in all calculated solutions without electron conductivity. 



However, numerical simulations by (jlgumenshchevl 120061 ) of non- rotating flows find ev- 
idence of convection. This convection may be physical. My model averages heat from all 
dissipation events over the fluid. Local reconnection events can lead to burst-type local 
heating that leads to buoyancy o f blobs. Also, magn etic buoyancy and diffusion play im- 
portant role in transfer processes (jlgumenshchevl 120061 ) . The correct inclusion of convection, 
magnetic buoyancy and diffusion is the subject of future studies. 
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=> stronger magnetic field, 
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(1c) outer magnetization <r„ 0d) outer magnetization o-„ 



Fig. 1. — Maximum accretion rate solution. Dependence of the accretion rate in units 
of Bondi rate on dimensionless parameters: characteristic length scale 7 (Fig. outer mag- 
netic helicity ^ (Fig. Ufa), outer matter magnetization (Fig. [Tb). Dependence (Fig. [TH) 
of sonic radius on outer magnetization a^. I take the reference model to have the following 
values of parameters: 7 = 1, = 1, ^ = 0.025. One parameter is varied to make one 
plot. Non-relativistic 1-T equation of state (dashed) verses relativistic 1-T equation of state 
(solid). 
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Isotropic outer turbulence turns Into anisotropic inner turbulence. 
Alfven point near outer boundary => no transport of angular momentum. 



0.00001 0.0001 0.001 0.01 0.1 1 

(2a) radius r/r B 




and inflow speed is smaller for non-relativistic EOS. 

0.00001 0.0001 0.001 o!oi OA 1 

(2b) radius r/r B 



Fig. 2. — Flow velocities, normalized to free-fall speed verses radius for maximum-rate 
solution: sound speed, inflow velocity, radial Alfven speed, 1-D perpendicular Alfven speed, 
turbulent velocity. Parameters — 1, 7 = 1, £00 = 0.025. Relativistic 1-T equation of state 
is on Figured, non-relativistic 1-T EOS is on Figure [2b. 
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Fig. 3. — Magnetization a verses dimensionless distance from the compact object r/rg is on 
Figure Dimensionless magnetic helicity £ verses dimensionless distance from the compact 
object t/tb is on Figure [3b. Both are for the maximum-rate solution with relativistic 
equation of state. 
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larger accretion rate M 
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Some deviations for relativistic EOS 
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Fig. 4. — Solution with angular momentum transport. Dependence of the accre- 
tion rate in units of Bondi rate on dimensionless parameters: characteristic length scale 7 
(Fig. Ek), outer magnetic helicity ^ (Fig. Hp), outer magnetization (Fig. Hp), and circu- 
larization radius r c ; r in units of r g (Fig. HH). I take the reference model to have the following 
parameters: 7 = 1, &00 = 1, ^cir = 10 3 r g , £00 = 0.025. Non-relativistic 1-T equation of state 
(dashed) verses relativistic 1-T equation of state (solid). 



-50- 




0.001 0.003 0.01 0.03 0.1 0.3 1 
(5a) radius r/r B 




0.001 0.003 0.01 0.03 0.1 0.3 1 
(5b) radius r/r B 



Fig. 5. — Flow velocities, normalized to free-fall speed verses radius for solution with 
angular momentum transport: sound speed, inflow velocity, radial Alfven speed, 1-D 
perpendicular Alfven speed, turbulent velocity. Parameters (Too = 1, 7 = 1, £ = 0.025, 
r cir = 10 3 r g . Relativistic 1-T EOS on Figured, non-relativistic 1-T EOS on Figure Eb. 
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Fig. 6. — Magnetization a verses dimensionless distance from the compact object rjr^ is on 
Figure [6k. Dimensionless magnetic helicity £ verses dimensionless distance from the compact 
object t/tb is on FigureEb. Both are for solution with angular momentum transport. 
Circularization radius is r r ; r = 10 3 r„. 
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Fig. 7. — Scheme of convection. Large magnetized blob is in perpendicular and radial 
pressure balance. Energy does not dissipate inside the blob. 



